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实 蝇 科 果 实 晶 属 昆虫 数字 图 像 目 动 识 别 系 统 
的 构建 和 测试 
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摘要 : 针对 双 翅 目 实 晶 科 果实 晶 属 昆虫 的 目 动 识 别 , 本 文 提 出 利用 起 及 中 脑 背 板 图 像 的 局 部 二 进 制 模式 (local 
binary pattern，LBP) 特征 ， 采 用 Adaboost 算法 ,设计 和 开发 “ 实 蝇 科 果 实 蝇 属 昆虫 数 字 图 像 目 动 识 别 系 统 ” 
(Automated Fruit fly Identification System-Bactrocera ，AFIS-B ) 。 该 系统 包括 图 像 采集 、 图 像 裁 前 、 预 处 理 、 特 征 提 
取 、 分 类 上 句 设 计 、 识 别 和 显示 , 共 7 个 模块 。 研 究 结 果 表 明 : LBP 特征 可 以 有 效 鉴 别 实 晶 科 果 实 蝇 属 昆虫 ; 在 对 实 
蝇 科 果实 蝇 属 8 个 种 的 测试 中 , 该 系统 表现 出 较 高 的 准确 性 和 稳定 性 , 平均 识别 率 可 达 80% 以 上 。 此 外 , 还 对 果 
实 蝇 属 昆虫 翅膀 及 中 脑 背 板 图 像 在 光照 不 均匀 、 姿 态 扭曲 、 样 本 受 损 及 样本 量 大 小 等 不 同 条 件 下 的 识别 率 进 行 了 
试验 测试 。 结 果 表 明 , 该 系统 对 测试 样本 的 光照 不 均 义 、 姿 态 扭曲 和 样本 受 损 都 表现 出 良好 的 鲁 棒 性 , 正确 识别 
率 与 训练 集 样本 各 个 种 数量 在 一 定 条 件 下 明显 正 相 关 , 与 训练 集 样 本 物种 总 量 负 相关 。 该 项 研究 为 实 蝇 科 有 害 昆 
虫 自动 识别 系统 的 构建 及 实际 应 用 提供 了 理论 、 方法 及 基础 数据 的 支撑 , 亦 可 为 其 他 昆虫 自动 识别 系统 的 研究 和 
构建 提供 有 益 借鉴 。 
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Abstract: Based on Local Binary Pattern (LBP) features of wing and scutum images and the improved 
Adaboost algorithm, we developed “ Automated Fruit Fly Identification System-Bactrocera, AFIS-B” for 
automatic identification of Bactrocera Macquart ( Diptera: Tephritidae ). The system consists of seven 
modules, which includes image acquisition, image cropping, image preprocessing, feature extraction, 
classifier design ，taxa identification and outcome display. The results showed that LBP features are 
effective to the automatic identification of fruit flies. The AFIS-B system has good accuracy and 
robustness by identifying 8 Bactrocera spp., and the average recognition rate ls more than 80%. We also 
did preliminary experiments under different conditions, such as inhomogeneous illumination, distorted 
posture, specimen partly damaged and different sample sizes. The results showed that the system has 
good robustness for the first three conditions, and the recognition rate usually positively relate to numbers 
of training sets for each species and negatively relate to the total species numbers. This research provides 
the theoratical , method and data foundation for the construction and practice of automated identification 
system of fruit fly, and it also gives a reference to the research and construction of other insects automated 
identification systems. 
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identification system 


长 久 以 来 , 生物 物种 的 鉴别 在 生物 学 领域 一 耳 
非常 重要 。 然 而 ,当今 世界 范围 内 传统 生物 分 类 学 
家 的 数量 却 在 不 断 下 降 。 如 何在 生物 分 类 学 家 人 数 
减少 的 同时 , 能 够 对 大 量 分 类 学 资源 进行 快捷 有 效 
的 信息 提取 ， 是 今天 我 们 面临 的 一 个 重要 挑战 
( Gaston and 0” Ne,，2004 ) 。 如 何 化 解 这 一 难题 和 
困境 ,可 能 会 一 定 程度 依赖 于 生物 物种 能 在 多 大 范 
围 内 实现 自动 化 的 识别 ( 陈 小 琳 等 , 2008 ) 。 现 代 计 
算 机 技术 的 高 速 发 展 ( 如 处 理性 能 的 飞速 提高 、 并 
行 计算 技术 的 成 熟 )、 模 式 识别 理论 的 日 益 完 善 、 
大 规模 数据 库 系统 的 成 功 应 用 、 网 络 和 数据 检索 技 
术 的 使 用 等 , 为 我 们 开发 方便 、 友 好 、 快 捷 的 物种 
分 类 工具 , 进行 大 规模 的 物种 自动 识别 提供 了 良好 
的 技术 平台 。 目 动 识别 系统 的 发 展 和 逐渐 普及 将 可 
能 让 分 类 学 家 从 繁杂 的 分 类 任务 中 解脱 出 来 , 专注 
于 新 物种 的 定义 和 发 现 等 , 并 为 进化 与 分 类 研究 、 
物种 多 样 性 研究 、 有 害 生物 的 监测 和 检疫 等 诸多 领 
域 提供 新 的 物种 识别 方法 。 

目前 , 研究 者 们 已 经 开发 出 了 一 些 基 于 图 像 或 
者 视频 的 物种 自动 识别 系统 ， 如 伦敦 自然 历史 博物 
馆 与 哥斯达黎加 大 学 ( 加利福尼亚 大 学 协会 ) 联合 
开发 的 基于 人 工 神 经 网 络 的 DAISY ( Digital 
Automated Identification System ) ( 0 ”Neill ef al., 
2000),， 中 国 农 业 大 学 沈 佐 锐 每 针对 棉铃 忠 
玉米 蜡 Pyrausta 
nubilalis (Hiibner ) 和 黑 哎 狂 蛙 Ectomocoris airox 
(St1) 等 设计 的 分 类 系统 (于 新 文 等 ,2003 ) 。 这 些 
基于 图 像 的 自动 识别 系统 一 般 都 是 针对 特定 类 群 或 
者 某 些 类 群 的 独 有 特征 而 设计 的 。 

实 晶 科 (Tephritidae) 隶 属于 双 功 目 ( Diptera ) 实 
蝇 总 科 (Tephritoidea ) 。 该 科 的 不 少 种 类 具有 十 分 
重要 的 经 济 意 义 , 其 中 果实 晶 属 Bactrocera Macquart 
包含 不 少 水 采 和 瓜 类 作物 的 危险 性 星 果 害 虫 , 如 梯 
小 实 晶 Bucirocera dorsalis ( Hendel ) ， 瓜 实 蝇 B. 
cucurbitae ( Coquillett ) 和 昆 干 兰 果 实 量 B. tryoni 
( Froggatt) 等 。 它 们 长 期 以 来 一 直 被 不 少 国 家 和 地 
区 列 和 人 严防 传人 的 检疫 性 有 害 生 物 名 单 中 。 但 在 世 
界 范 围 内 全 面 而 系统 地 开展 实 晶 科 昆虫 日 劲 识别 系 
统 的 研究 则 尚未 进行 。 基 于 此 背景 , 我 们 研制 了 针 
对 果实 蝇 属 昆虫 图 像 的 目 动 识别 系统 (Automated 


Fruit Fly ldentification System-Bactrocera Macquart, 


Helicoverpa armigera Hiibner, 


AFIS-B), 并 对 该 类 群 图像 的 目 动 识别 进行 了 较为 
深入 的 试验 研究 和 讨论 。AFIS-B 设计 的 初衷 是 用 
来 研究 和 自动 识别 果实 蝇 属 物种 , 但 将 来 则 有 可 能 
用 于 实 晶 科 其 他 类 和 群 的 目 动 识别 , 并 有 可 能 与 昆虫 
的 其 他 目 动 识别 系统 整合 , 从 而 为 生物 的 目 动 分 类 
与 识别 等 提供 更 多 便利 。 


1 实 蝇 科 果 实 蝇 属 昆 虫 图 像 自 动 鉴 别 
系统 (AFIS-B ) 的 构建 


1.1 开发 工具 与 平台 

便 件 平台 : 明美 ME21 数字 图 像 采 集 系统 一 套 
( 奥 林 巴 斯 CX2 显微镜 , 500 万 像素 数码 摄像 系统 )、 
CANON PowerShotA720 IS 数码 相机 一 台 、PC 机 一 合 
(Intel core 2 Duo T6600, 2.2GHz, 2.0G 内 存 )。 

软件 平台 : Microsoft SQLServer2005 ，Matlab7.0。 

开发 工具 : VC ++6.0。 
1.2 系统 结构 设计 

AFIS-B 的 总 体 设 计 思 想 是 ,通过 对 果实 蝇 属 
昆虫 各 类 具有 分 类 学 意义 图 像 的 研究 , 选取 适 于 计 
算 机 进行 目 动 识别 的 特征 , 利用 训练 集 样本 进行 图 
像 特征 提取 和 分 类 带 设 计 , 用 测试 集 样本 测试 其 目 
动 识别 的 正确 率 。 

系统 共 包 含 7 个 模块 : 图 像 采 集 、 图 像 区 域 裁 
前 、 预 处 理 、 特 征 提取 、 分 类 副 设 计 、 识别 和 显示 
模块 , 各 模块 之 间 的 关系 见 图 1。 
1.3 系统 模块 功能 介绍 
1.3.1 图 像 采 集 : 果实 蝇 的 翅 是 人 工 鉴别 时 使 用 
的 主要 特征 之 一 , 其 图 像 有 明显 的 怒 脉 和 翅 斑 , 日 
基本 处 于 同一 平面 ,因此 非常 适合 用 来 做 自动 识 
别 ; 中 胸 背 板 与 翅膀 相 比 , 立体 性 稍 强 , 但 也 相对 
较 遍 平 ,， 并 具有 规律 的 条 带 和 斑纹 , 因而 也 是 人 工 
鉴别 时 常 使 用 的 重要 特征 之 一 。 我 们 使 用 数字 图 像 
采集 系统 ,对 标本 的 这 膀 及 背 板 图 像 进行 了 重点 
采集 。 

采集 时 显微镜 连接 数码 摄像 头 , 摄像 头 则 用 数 
据 线 与 计算 机 相连 , 将 符合 规格 的 标本 按 正 确 的 角 
度 放 置 在 显微镜 下 , 并 聚焦 于 采集 部 位 的 中 心 。 在 
计算 机 上 获得 放大 的 数字 图 像 后 , 将 图 像 保 存 成 数 
字 图 像样 本 。 男 外 ,也 可 直接 用 数码 相机 对 准 显 微 
镜 目镜 进行 数字 图 像 的 简易 采集 。 
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a 图 像 采 集 b 图 人 裁 前 小 c 图 像 预 处 理 
从 数码 相机 或 者 图 像 | 图 像 尺 寸 的 统一 以 及 彩色 图 像 
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图 1 实 蝇 科 果实 蝇 属 昆虫 图 像 目 动 鉴别 系统 (AFIS-B) 各 模块 关系 图 
Fig. 1 The flow chart of seven modules of AFIS-B 
a: 图 像 采 集 Image acquisition (从 数码 相机 或 者 图 像 文 件 等 获得 Acquire from digital camera or image file) ; b: 图 像 裁 前 Image cropping (得 到 需 
要 识别 部 分 Obtain recognized part) ; c: 图像 预 处 理 Image preprocessing ( 图像 尺 寸 的 统一 以 及 彩色 图 像 灰 度 化 、 归 一 化 Make image size uniform 
and transform colored images into gray levels ) ; d: 特征 提取 Feature extraction (从 图 像 中 提取 适 于 分 类 的 特征 Extract classification features from 
image ) ; e: 分 类 器 设计 Classifier design (调整 分 类 器 参数 Adjust classifier parameter) ; f: 识别 模块 Identifying module (识别 输入 图 像 Identify the 
input image) ; g: 显示 模块 Displaying module (显示 识别 结果 Display the recognized result ) . 


1.3.2 ”图像 裁剪 与 预 处 理 : 利用 本 系统 的 “图 像 区 ”将 和 矩形 框 中 的 图 像 裁 前 并 保存 , 并 用 预 处 理 模 块 先 
域 裁 冀 模块 ”中 的 图 像 裁 前 工具 进行 翅 及 背 板 部 分 。 把 裁 冯 出 来 的 图 像 转换 成 灰 度 图 , 再 进行 图 像 归 一 
的 裁 前 。 裁 艾 工 具 的 原理 是 先 把 图 像 经 过 旋转 调整 。” 化 处 理 , 以 提高 系统 的 处 理 速 度 和 识别 率 , 其 主要 
到 统一 的 姿态 , 再 利用 鼠标 在 图 像 上 设 定 矩 形 框 ， ”流程 见 图 2。 


B 
图 2 翅 图 像 区 域 裁剪 及 预 处 理 


Fig. 2 The cropping and preprocessing of wing image 
A: 原始 图 像 翅膀 区 域 选 择 Select the original area of wing image; B: 裁剪 出 的 翅膀 图 像 The cropped wing; C: 灰 度 化 后 的 翅膀 图 像 The wing 


transformed into gray levels. 
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1.3.3 图 像 特征 提取 : 为 了 实现 多 类 鉴别 技术 ， 
系统 采用 了 正 负 样 本 的 思想 ( 叶 剑 华 和 刘 正 光 ， 
2008 ) ,通过 不 同 的 LBP 算 子 和 滑动 窗 技 术 , 间接 
进行 样本 特征 的 提取 。 首 先 通过 滑动 窗 获 得 两 幅 样 
本 对 应 区 域 , 然后 计算 对 应 区 域 每 个 像素 的 局 部 二 





进 制 模式 (local binary pattern，LBP) 值 (Ojala et al.， 
2002), 最 后 计算 对 应 区 域 LBP 直方 图 相似 值 (x 
距离 ), 并 将 此 值 作为 一 个 特征 。 通 过 滑动 窗口 和 
改变 LBP 算 子 来 获得 高 维 的 特征 向 量 , 具体 见 
图 3。 





Dn 
S(H,, H)= SB (HH AH tH) 
了 


图 3 通过 滑动 窗 和 LBP 算 子 来 获取 一 个 特征 示例 
Fig. 3 Example of acquiring feature by using slide window and LBP operator 
a: 滑动 窗 A slide window; b: 区 域 LBP 直方 图 LBP histograms of corresponding windows; c: 两 个 直方 图 的 卡 方 距 离 The x? distance of two 


histograms. 


1.3.3.1 LBP 算 子 : LBP 算 子 用 于 描述 图 像 中 每 
个 像素 值 与 其 邻 域 像素 值 的 大 小 关系 , 在 灰 度 图 像 
中 , 像素 值 为 像素 的 灰 度 。 基 本 的 LBP 算 子 可 用 图 
4 进行 描述 。 将 3 x3 邻 域 周围 8 个 像素 与 中 心 像 
素 值 进行 比较 并 二 值 化 , 按 顺 时 针 方 向 排列 得 到 一 
个 8 位 二 进 制 数 , 再 将 其 转化 为 十 进 制 数 , 将 该 值 
作为 中 心 像素 的 编码 , 该 编码 就 是 LBP 算 子 对 中 心 
像素 作用 后 的 结 
二 进 制 Binary 
11010001 


十 进 制 Decimal 
209 











回 马 听 司 二 星 前 阿 
Threshold 
| 


图 4 基本 的 LBP 算 子 示例 
Fig. 4 Example of simple LBP operator 


LBP 可 用 不 同 尺 寸 的 邻 域 来 表示 。 采 用 圆 形 邻 
域 结合 双 线 性 插值 ， 邻 域 的 半径 和 像素 数目 可 任意 
选取 。LBP 用 (P, R ) 表 示 半 径 为 R 的 贺 上 的 P 个 
采样 点 。 图 5 是 圆 形 的 (8, 2) 邻 域 ， 而 基本 的 LBP 
算 子 则 是 (8，1 ) 邻 域 。 





图 5 圆 形 (8, 2) 上 的 LBP 采样 点 
Fig. 5 LBP sample point on circular (8, 2) neighborhood 
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图 像 LBP 编码 的 直方 图 H; 可 定义 如 下 : 
H, = >》 Tf(x,y) = 村 ,zz=0,…7 -1 (1) 


其 中 , n 是 LBP 算 子 产生 的 不 同 编码 数量 。 如 尸 = 
8, 则 =256。 及 x,y) 是 坐标 x , y 处 的 编码 。 而 : 
1, 4 为 真 
(4 = {0 4 pm (2) 
1.3.3.2 区 域 LBP 下 方 图 : 式 (1) 的 直方 图 包含 
了 整个 图 像 上 的 微 结构 分 布 ,， 如 边 、 点 和 平滑 区 
域 。 为 了 更 好 地 表示 图 像 , 应 该 保留 图 像 的 空间 结 
构 信 息 。 为 此 , 在 图 像 上 放置 一 个 窗口 , 并 进行 平 
移 ( 图 3: a) ,可 得 到 区 域 Ro , Ri ，… , R,, ， 而 区 
域 LBP 直方 图 (regional LBP histogram ,RLBPH ) (图 
3: b) 可 定义 为 : 
全 一 >》 Tf(x,y) = jj} (x,y) E R,| ， 
7] = 0,…,n 1 =0,…,m—1 (3) 
而 两 个 直方 图 的 卡 方 距离 (图 3: c) 则 定义 为 : 


S(Hi,E) = (Hi Hi)/ (Hi + Hi) (4) 


其 中 : n 是 LBP 算 子 产生 的 不 同 编码 数量 。 如 P = 
8 ， 则 n=256。H; , H, 是 两 个 ( 区 域 ) 直方 图 ， Hi;、 
五; 分别 是 五 和 五 第 i 个 编码 出 现 的 次 数 。 此 距离 
(区 域 直方 图 相似 度 值 ) 即 可 作为 样本 的 一 个 特 
征 值 。 

1.3.4 分 类 器 设计 : 通过 滑动 窗 的 不 断 平 移 ，1. 
3.3.2 中 的 每 个 样本 的 区 域 耳 方 图 相似 度 值 (特征 
值 ) 可 达 上 千 个 。 如 此 众多 的 特征 包含 大 量 的 元 余 
言 息 , 因此 ,系统 分 类 器 设计 采用 了 将 弱 分 类 器 融 
合成 强 分 类 器 的 方法 。 训 练 算法 是 AdaBoost 
( Freund, 1995; Freund and Schapire，1999 ) ， 即 从 
高 维 的 特征 问 量 中 提取 出 分 类 准确 率 最 高 的 厂 干 的 
子 回 量 来 设计 奉 干 弱 分 类 天 ,， 并 赋 子 每 个 弱 分 类 天 
不 同 的 权重 (权重 与 该 弱 分 类 带 的 分 类 正确 率 正 相 
关 ) , 通过 算数 和 的 方式 将 弱 分 类 融融 合 起 来 ， 最 
后 的 分 类 绪 采 是 所 有 弱 分 类 需 输 出 结果 的 加 权 算 
术 和 。 

1.3.4.1 AdaBoost: AdaBoost 是 组 合 分 类 大 最 常用 
的 算法 , 它 逐 步 将 弱 分 类 副 线 性 组 合成 强 分 类 兹 。 
当 用 于 特征 选择 时 ，AdaBoost 基于 选择 的 某 个 特征 
构造 弱 分 类 璐 ,最 终 的 强 分 类 需 是 这 些 弱 分 类 闪 的 
线性 组 合 。 

1.3.4.2 多 类 鉴别 技术 : AdaBoost 本 质 上 是 通过 
学 习 得 到 二 类 分 类 器 。 为 了 将 多 类 问题 转化 为 二 类 
问题 , 本 文 借鉴 了 类 内 差 和 类 间 差 的 思想 。 这 里 的 


类 内 差 由 同 种 的 两 个 不 同 图 像 的 多 个 RLBPH 相似 
度 构成 ， 类 间 差 由 不 同 种 的 两 个 图 像 的 多 个 
RLBPH 相似 度 构成 ， 即 可 得 到 N 个 区 域 直方 图 相 
似 度 值 作为 特征 值 序列 ，N 的 选取 取决 于 选择 的 滑 
动 窗 的 大 小 和 步 长 。 

1.3.4.3 训练 测试 过 程 : 对 于 给 定 的 训练 集 , 类 
内 差 ( 正 例 ) 样 本 和 类 间 差 (反例 ) 样 本 数目 通常 相 
差 很 大 。 例 如 训练 集中 有 4 个 种 ,每 种 个 样本 。 
由 作 差 法 形成 的 正 例 和 反例 样本 集 数目 分 别 为 : N' 
= CC; = gk(k—-1)/2, 和 N = CCiC; = kgq(g 
- DA 。 它 们 的 比值 为 : -= “站 一。 假定 4 


=100,k=10, 则 有 N 和 AN 分 别 为 49$ 000 和 
4 500。 反 例 样 本 集 数 目 和 正 例 样 本 集 数 目 比 值 
为 110。 

本 文采 用 AdaBoost 算法 分 别 从 灰 度 图 像 的 众 
多 RLBPH 相似 性 特征 中 选择 最 有 利于 分 类 的 少量 
特征 ,构造 组 合 分 类 项 。 该 算法 的 流程 如 下 : 

输入 : 整个 训练 集 { (xn ,7 ) (xn yy) |， 
为 样本 , y;，e 1 -1,1| 为 标识 (1 表示 正 例 样本 ， 
-1 表示 反例 样本 )。 

初始 化 : 设 训练 集中 所 有 的 正 例 样本 数目 为 
NM ， 负 例 样 本 数目 为 为 V 。 





Begim 
1) 初 始 化 训练 样本 的 权重 ， 即 : 
1 
yi = 1 
Wo,; 一 " 1 ) | (5) 
CN 


Fori =1,.…,(N+N ) 。 

2) 在 当前 的 训练 集 上 针对 某 个 特征 S,(x) , 通 
过 简单 取 国 值 0 得 到 弱 分 类 帮 h(x) : h,(x) = 
sign(S,(%) — 0) (6) 

弱 分 类 器 的 选择 以 使 本 次 分 类 错误 率 最 低 的 弱 
分 类 融 作 为 弱 分 类 囊 。ow 为 h,(x) 的 权重 , 其 值 为 : 





1, ,1l-e, 
wu = In( 2) (7) 
其 中 : s, 为 弱 分 类 央 的 分 类 错误 率 。 
(N++N-) 
2; = 之 w, ;8 六 (Xi ) 天 y;) (8) 
3 ) 更改 样本 的 权重 
Ww, ; 本 让 h(%;) = yy 
Wir1 一 x 一 
Ls [ew if h(x) 天 
1 Ca 
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其 中 : Z, 是 归 一 化 因子 , 以 使 w, 仍然 为 某 种 分 布 ， 
本 3 W, ; (10) 
4) 组 合 h, 和 以 前 学 习 的 强 分 类 带 Hi 形成 强 

分 类 融 有， 


1 = sien(H, (x) +h,) = 


sign( Da «hls)) (11) 


上 述 AdaBoost 学 习 算 法 , 可 获得 区 分 类 内 差 和 
类 间 差 的 强 分 类 髓 。 因 为 AdaBoost 学 习 到 的 弱 分 
类 器 对 应 RLBPH 相似 性 特征 ， 故 本 文 将 w 作为 
RLBPH 相似 度 特征 的 权重 , 采用 式 (4) 计 算 两 幅 图 


像 的 相似 度 特征 。 则 这 两 个 局 部 图 像 的 匹配 得 
分 为 ; 
Eh) 1 
给】 三 1 2 
fAR) Ne I (12) 


这 两 幅 ee 7 人 为 : 
(4) - = Yah) (13) 
T 为 强 分 类 器 含有 的 弱 分 类 器 数量 。 该 方法 称 
为 AdaRLBPH 。 
1.4 系统 界面 及 主要 功能 
1.4.1 系统 主 界面 (图 6): 
系统 主 界面 主要 包括 供用 户 选 择 使 用 的 识别 类 
型 : 计 膀 识别 、 背 板 识 别 、 腹 部 识别 、 交 互 式 检索 ， 


F(x) = 


或 者 是 选择 查询 寄主 信息 等 。 用 户 可 以 通过 图 标 或 
左 侧 的 栏目 名 称 进 入 下 一 级 界面 。 目 前 只 提供 翅膀 
识别 、 背 板 识 别 两 个 功能 模块 。 

1.4.2 翅膀 图 像 识别 子 系统 界面 (图 7): 该 界面 
的 主要 功能 是 训练 分 类 融 、 测 试 分 类 天 分 类 结果 、 
用 训练 出 来 的 分 类 需 测试 识别 率 和 输入 用 户 指定 的 
符 识 别 翅膀 图 像 ， 以 及 转 到 其 他 识别 子 系统 或 
查询 。 

1.4.3 翅膀 图 像 识 别 结果 显示 模块 (图 8) : 该 模 
块根 据 用 户 选 定 的 待 识 别 翅 膀 图 像 , 给 出 识别 结 
果 , 并 列举 出 相似 度 最 高 的 前 两 个 物种 图 像 。 符 
系统 进一步 完善 后 , 将 同时 显示 所 识别 物种 的 其 
他 形态 特征 、 地 理 分 布 、 生 物 学 和 经 济 意义 等 
信息 。 

1.4.4 中 胸 背 板 图 像 识 别 子 系统 界面 (图 9): 该 
界面 的 主要 功能 是 训练 分 类 右 、 测 试 分 类 带 分 类 结 
果 、 用 训练 出 来 的 分 类 天 测试 识别 率 和 输入 用 户 指 
定 的 待 识别 中 胸 背 板 图 像 ， 以 及 转 到 其 他 识别 子 系 
统 或 查询 。 

1.4.5 中 胸 背 板 图 像 识别 结 采 显示 模块 (图 10): 
该 模块 根据 用 户 点 选 的 待 识别 背 板 图 像 , 给 出 识别 
结果 , 并 列举 出 相似 度 最 高 的 前 两 个 物种 图 像 。 待 
系统 进一步 完善 后 , 将 同时 显示 所 识别 物种 的 其 他 
形态 特征 、 地 理 分 布 、 生 物 学 和 经 济 意义 等 信息 。 
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图 6 实 蝇 科 果实 蝇 属 昆虫 数字 图 像 自 动 识别 系统 (AFIS-B) 主 界面 
Fig. 6 The main interface of Automated Fruit Fly Identification System-Bactrocera (AFIS-B ) 
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实 蝇 科 果 实 蝇 属 昆虫 数字 图 像 自 动 识 别 系统 (AFIS-B) 翅 膀 识 别 子 系统 界面 


Fig. 7 The interface of wing image identification subsystem of AFIS-B 
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Distribution Unknown 
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图 8 翅膀 图 像 识别 结果 显示 模块 
Fig. 8 The recognization result module of wing image identification subsystem of AFIS-B 
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图 9 实 蝇 科 果 实 蝇 属 昆虫 数字 图 像 自 动 识 别 系 统 (AFIS-B) 中 胸 背 板 识别 子 系统 界面 
Fig. 9 The interface of thorax image identification subsystem of AFIS-B 
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图 10 中 胸 背 板 图 像 识 别 结 采 显 示 模 块 
Fig. 10 The recognization result module of thorax image identification subsystem of AFIS-B 


2 识别 试验 测试 结果 


本 文 针对 果实 蝇 属 昆虫 标本 中 胸 背 板 或 翅膀 图 
像 在 光照 不 均 义 、 姿 态 扭曲 、 样 本 受 损 及 样本 量 大 
小 等 不 同 条 件 下 的 识别 率 进 行 了 试验 测试 。 所 有 试 
验 用 到 的 物种 的 测试 集 样本 均 不 与 训练 集 样本 重 
合 。 测 试 集 测试 的 时 候 , 判断 某 幅 待 测 图 像 属于 茶 


个 物种 的 依据 为 : 该 图 像 匹 配 到 某 物 种 某 个 图 像 的 
得 分 最 高 。 分 类 句 训 练 和 测试 时 的 特征 选择 所 采用 
的 LBP 滑动 窗口 大 小 为 15 pixel x 15 pixel, 水 平 垂 
直方 同 的 滑动 步 长 均 为 1 pixel。 识 别 率 是 指 物种 测 
试 集 正确 识别 出 的 样本 量 与 参加 测试 的 物种 样本 量 
的 比值 ; 总 识别 率 是 指 几 个 物种 测试 集 正确 识别 出 
的 样本 量 与 参加 测试 的 总 样本 量 的 比值 。 
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2.1 试验 材料 

采集 实 蝇 科 果实 晶 属 昆虫 14 个 物种 的 标本 图 
像 , 与 本 文 研究 有 关 的 物种 图 像 采集 信息 见 表 1。 
并 选取 其 中 8 个 物种 的 翅 或 中 胸 背 板 图 像 作为 研究 
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对 和 象 (物种 名 、 相 关 样 本 的 数量 均 已 在 表 1 加 粗 标 
出 ) ， 对 每 个 物种 3 ~ 10 个 以 上 个 体 的 图 像 进 行 训 
练 和 测试 。 所 有 训练 和 测试 集中 的 样本 均 是 基于 分 


表 1 果实 蝇 属 昆虫 图 像 采 集 信息 表 
Table 1] Image information of Bactrocera Macquart 


充 膀 Wing 中 胸 背 板 Scutum 
下 物种 名 年 本 本 。 完好 图 像 量 。 受 损 图 像 量 ”姿态 扭曲 图 像 ”样本 量 ( 幅 ) 。 完好 图 像 ns se 
Index Species name Pa ( 幅 ) ( 幅 ) 量 ( 幅 ) Number of 量 ( 幅 ) 0 
Number of Inhomogeneous 
ee Intact samples Damaged samples Distorted samples samples Intact samples OE 
1 B. (B. ) carambolae 9 6 1 1 4 4 0 
2 B. (B. ) dorsalis 25 12 12 6 83 33 S0 
3 B. (B. ) latifrons 4 0 4 3 1 1 0 
4 B. (B. ) rubigina 20 13 4 3 77 14 63 
5 B. (P. ) garciniae 5 2 之 3 4 4 0 
6 B. (T. ) minax 6 5 1 lL 6 6 0 
7 B. (T. ) tsuneonis 13 5 2 6 4 4 0 
8 B. (2Z. ) abbreviata 10 6 2 1 3 3 0 
9 B. (2Z. ) caudata 46 6 7 18 42 12 30 
10 B. (2Z.) cucurbitae 101 17 78 67 78 S3 25 
11 B. (Z.) diaphora 8 4 4 5 5 0 
12 B. (Z. ) scutellaris 3 3 0 0 4 4 0 
13 B. (2Z. ) scutellata 18 11 6 3 1S 1S 0 
14 B. (Z.) tau 103 21 8 31 21 21 0 


表 中 样本 物种 名 称 的 属 名 和 亚 属 名 均 进 行 了 如 下 缩写 The genus and subgenus names are abbreviated as below: B. = Bactrocera, P. = Paratridacus,， 
也 = Tetradacus, Z. = Zeugodacus， 加 粗 行 表 示 本 文 所 做 的 相关 识别 率 测 试 实验 均 以 这 些 样本 为 测试 对 象 。The species in bold are the samples 


used in testing experiments. 


2.2 光照 不 均匀 对 识别 率 的 影响 

光照 不 均匀 图 像 是 指 采集 图 像 时 , 样本 整体 或 
者 部 分 受 光 不 均匀 ,致使 获取 到 的 数字 图 像 的 不 同 
部 分 亮度 不 一 ， 此 试验 以 中 胸 背 板 图 像样 本 作为 研 
究 对 象 (图 11) 。 分 类 需 训 练 和 测试 结果 见 表 2。 

表 2 列 出 了 每 种 测试 集 各 种 样本 (正常 、 光 照 
不 均匀 、 全 部 ) 的 识别 率 ( 第 1~4 行 ,7~9 列 ) 和 总 
识别 率 ( 第 5 行 ,7~9 列 )。 第 5 行 的 第 3 ~6 列 则 
是 每 列 图 像 数 量 之 和 。 研 究 结果 表明 , 在 训练 集 物 
种 数 为 4 的 情况 下 , 每 种 正常 样本 的 识别 率 均 高 于 
光照 不 均匀 样本 的 识别 率 , 但 是 差异 不 是 很 大 。 由 
此 验证 了 LBP 特征 对 光照 不 均匀 所 造成 的 影响 不 
2.3 姿态 扭曲 对 识别 率 的 影响 

姿态 扭曲 图 像 是 指 图 像 整 体 或 者 部 分 产生 形 
变 , 与 图 像 的 正常 形状 产生 了 偏差 后 采集 到 的 图 


像 , 此 试验 用 翅膀 图 像样 本 作为 研究 对 象 (图 12)。 
分 类 带 训练 和 测试 结果 见 表 3。 


pr 





A B 
图 11 物种 正常 图 像样 本 (A) 与 光照 不 均匀 图 像 
样本 (B) 的 比较 
Fig. 11 Comparison of normal samples (A) and 


inhomogeneous illumination samples (B) 
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表 2 光照 不 均匀 试验 测试 结果 


Table 2 Result of testing experiment of Bactrocera Macquart under inhomogeneous ilumination 


训练 集 样本 量 ( 幅 ) 


Training set samples 


测试 集 样本 量 ( 幅 ) 


Testing set samples 


识别 率 (% ) 


Rate of identification 





品 人 
Nise 正常 图 像 ”光照 不 均匀 图 像 ”正常 图 像 ”光照 不 均匀 图 像 ”正常 图 像 ” ”光照 不 均匀 图 像 ”全 部 图 像 
Normal Inhomogeneous Normal Inhomogeneous Normal Inhomogeneous Total 
samples illumination samples samples illumination samples samples illumination samples samples 
1 B. (B. ) dorsalis 9 8 42 24 97.6 87.5 93.9 
2 B. (B. ) rubigina 4 9 54 10 85.2 80.0 84.4 
3 B. (Z. ) caudata 4 7 23 8 95.7 87.5 93.5 
4 B. (Z. ) cucurbitae 12 6 19 41 100 87.8 91.7 
总 计 Total 29 30 138 83 92.8 86.8 90.5 
图 12 物种 正常 图 像样 本 (A) 与 姿态 扭曲 图 像样 本 (B) 的 比较 
Fig. 12 Comparison of normal samples (A) and distorted posture samples (B) 
表 3 姿态 扭曲 试验 测试 结果 
Table 3 Result of testing experiment of Bactrocera Macquart with distorted posture samples 
训练 集 样本 量 ( 幅 ) 测试 集 样 本 量 ( 幅 ) 识别 率 (% ) 
Training set samples Testing set samples Rate of identification 
序号 物种 名 
和 正常 图 像 扭曲 图 像 ” ”正常 图 像 ” ”扭曲 图 像 正常 图 像 扭曲 图 像 ” ”全 部 图 像 
Normal Distorted Normal Distorted Normal Distorted Total 
samples posture samples samples posture samples samples posture samples samples 
1 B. (Z. ) caudata 之 4 4 14 50.0 78.6 72.2 
2 B. (Z. ) cucurbitae 5 9 12 58 91.7 93.1 92.9 
3 B. (Z.) tau 7 5 14 26 71.4 37.7 62.5 
总 计 Total 14 18 30 98 76.7 81.6 80.5 


表 3 列 出 了 每 种 测试 集 各 种 样本 (正常 、 扭 曲 、 
全 部 ) 的 识别 率 ( 第 1~3 行 ,7~9 列 ) 和 总 识别 率 
(第 4 行 ,7~9 列 )。 第 4 行 的 第 3 ~6 列 则 是 每 列 
图 像 数 量 之 和 。 结 果 表明 , 在 训练 集 物种 数 为 3 的 
情况 下 , 每 种 正常 样本 的 识别 率 与 姿态 扭曲 样本 的 
识别 率 的 差异 均 不 是 很 明显 LB. (2Z. ) caudata 除 
外 , 因为 当 正 稼 图 像 训练 集 样本 量 较 少 时 , 这 是 一 
种 概率 现象 无 法 避免 ], 这 可 能 是 因为 LBP 特征 具 
有 旋转 不 变性 的 特点 。 同 时 ，Adaboost 算法 在 选择 
特征 生成 弱 分 类 带 的 时 候 , 避 开 了 选择 姿态 扭曲 部 


分 , 而 选择 了 分 类 准确 率 最 高 的 特征 。 因 此 ,影响 
识别 率 的 主要 因素 集中 在 了 训练 集 样本 量 的 大 小 上 
面 。 由 此 验证 了 此 系统 对 姿态 扭曲 图 像 的 影响 不 
敏感 。 
2.4 标本 受 损 对 识别 率 的 影响 

标本 受 损 图 像 是 指标 本 受到 轻微 损伤 ( 指 受 损 
面积 不 超过 总 面积 的 1/3, 且 保 留 主要 识别 特征 ) 
后 采集 到 的 图 像 ， 此 试验 用 翅膀 图 像 作 为 样本 (图 
13)。 分 类 右 训 练 和 测试 结果 见 表 4。 
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图 13 物种 正常 图 像样 本 ( A) 与 受 损 图 像样 本 (B) 的 比较 
Fig. 13 Comparison of normal samples ( A) and damaged samples (B) 
表 4 标本 受 损 试验 测试 结果 
Table 4 Result of testing experiment of Bactrocera Macquart with slightly damaged specimen 
训练 集 样本 量 ( 幅 ) 测试 集 样本 量 ( 幅 ) 识别 率 (% ) 
Training set samples Testing set samples Rate of identification 
序号 物种 名 
index 。 Species name 正常 图 像 。”” 受 损 图 像 ”正常 图 像 ”” 受 损 图 像 ”正常 图 像 受 损 图 像 ” ”全 部 图 像 
Normal Damaged Normal Damaged Normal Damaged Total 
samples samples samples samples samples samples samples 
1 B. (Z. ) caudata 必 3 + 4 100 100 100 
2 B. (Z. ) cucurbitae 5 10 12 68 91.7 97.1 96.3 
3 B. (Z.) tau 党 3 14 3 92.9 100 94.7 
总 计 Total 14 16 30 77 93.3 97.4 96.3 


表 4 列 出 了 每 种 测试 集 各 种 样本 (正常 、 受 损 、 
全 部 ) 的 识别 率 ( 第 1~3 行 ,7~9 列 ) 和 总 识别 率 
(第 4 行 ,7~9 列 )。 第 4 行 的 第 3 ~6 列 则 是 每 列 
图 像 数 量 之 和 。 研 究 络 采 表明 , 在 训练 集 物种 数 为 
3 的 情况 下 , 每 种 正常 样本 的 识别 率 均 与 标本 受 损 
样本 的 识别 率 的 差异 不 是 很 明显 , 这 是 因为 
Adaboost 算法 在 选择 特征 生成 弱 分 类 天 的 时 候 ， 避 
开 了 选择 标本 受 损 部 分 , 而 选择 了 分 类 准确 率 最 高 


的 特征 。 由 此 验证 了 此 系统 对 标本 受 损 图 像 的 影响 
不 敏感 。 
2.5 每 种 训练 集 样本 数量 对 识别 率 的 影响 

本 试验 测试 每 种 训练 集 样本 数 对 识别 率 的 影 
啊 , 通过 寻找 识别 率 与 每 种 训练 集 样本 数 的 依赖 关 
系 , 确定 训练 集 最 佳 的 样本 个 数 。 此 试验 用 这 膀 图 
像 作为 样本 , 分 类 融 训 练 和 测试 结果 见 表 5。 


表 5 训练 集 每 种 样本 数量 试验 测试 结果 


TableS Result of testing experiment of Bactrocera Macquart with different training sample size for each Species 


每 种 识别 率 Identification rate of each species (%) 


训练 集 样本 量 ( 幅 ) 


Training set samples 


B. (B. ) dorsalis 


00) (7 
50 57.1 
3 66.7 85.7 
4 83.3 71.4 
5 50 85.7 
6 33.3 100 
7 66.7 85.7 
8 66.7 85.7 


* 每 种 测试 集 样本 数量 ( 幅 ) The testing set samples for each species. 


B. (B.) ruwbigina B. (Z.) cucurbitae 


i 总 识别 率 (% ) 
Total identification rate 
(10)“ (12) 
20 25 34.3 
30 58.3 S71 
80 66.7 74.3 
100 100 88.6 
80 91.7 80 
10 100 91.4 
90 100 88.6 


2 期 


表 5 列 出 了 训练 集 每 种 样本 数量 从 2 ~ 8 的 不 
同情 况 下 , 每 种 测试 集 样本 的 识别 率 ( 第 4 ~ 10 行 ， 
2 ~5 列 ) 和 总 识别 率 ( 第 6 列 ) 。 总 识别 率 是 指 每 行 
所 有 物种 图 像 的 测试 集 样本 的 总 识别 率 。 研 究 结 果 
表明 , 在 训练 集 物种 数量 为 4 的 情况 下 ， 随 着 每 种 
训练 集 样本 数量 的 增加 , 识别 率 呈现 明显 的 上 升 趋 
擅 ; 每 种 训练 集 样本 数量 为 7 的 情况 下 ,识别 率 最 


表 6 
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高 , 并 达到 90% 以 上 。 当 每 种 训练 集 样本 数量 超过 
7, 识别 率 不 再 明显 上 升 , 反而 有 所 降低 , 这 与 分 类 
人 船 过 学 习 有 天 。 
2.6 ”物种 数量 对 识别 率 的 影响 

本 试验 测试 建 库 物种 总 数量 对 识别 率 的 影响 ， 
探讨 系统 用 于 大 量 物种 识别 的 可 能 性 ,此 试验 用 翅 
膀 图 像 作为 样本 。 分 类 天 训练 和 测试 结果 见 表 6。 


物种 数量 试验 训练 测试 结果 


Table 6 Result of testing experiment of Bactrocera Macquart with different total species numbers 


训练 集 样本 测试 集 样本 


序号 物种 名 量 ( 幅 ) 量 ( 幅 ) 
Index Species name Training Testing 
set samples set samples “ spp 
1 B. (2Z.) tau 7 14 100 
2 B. (2Z.) cucurbitae 5 12 91.7 
3 BB. (B.) ruwbigina 5 8 — 
4 B. (B. ) dorsalis 4 8 — 
5 B. (Z. ) scutellata 4 7 一 
6 B. (2Z.) caudata 3 3 — 
7 B. (B.) carambolae 3 3 — 
8 B. (2Z.) abbreviata 3 3 — 
图 像 总 量 Total samples 34 58 26 
总 识别 率 (% ) 96 2 


Total identification rate 


表 6 列 出 了 物种 数量 在 2 ~8 的 不 同情 况 下 ， 
每 种 测试 集 样本 的 识别 率 (第 1~8 行 ,5~11 列 ) 
和 总 识别 率 ( 第 10 行 )。 图 像 总 量 在 第 3 ~4 列 是 
指 每 列 图 像 数 量 之 和 , 在 第 5 ~ 11 列 指 每 列 受 测 物 
种 测试 集 样本 图 像 数 量 之 和 ,总 识别 率 是 指 每 列 所 
有 受 测 物种 图 像 的 测试 集 样本 的 总 识别 率 。 研 究 结 
采 表 明 , 随 者 训练 集 物种 数 的 增加 , 识别 率 呈 现下 
降 趋 势 ， 由 于 物种 序号 为 4 及 以 后 的 物种 的 训练 集 
样本 量 偶 少 , 导致 识别 率 下 降 更 加 迅速 。 该 结论 给 
我 们 进行 海量 物种 识别 提供 了 依据 , 即 必 须要 使 每 
种 训练 集 样本 达到 合适 的 数量 。 

上 述 试验 结果 表明 : 光照 不 均 义 、 姿 态 扭曲 及 
样本 局 部 受 损 等 对 最 后 的 识别 结果 影响 不 是 很 明 
显 , 这 是 因为 LBP 特征 具有 光照 鲁 棒 性 、 旋 转 不 变 
形 , 而 且 分 类 副 进 行 特征 选取 的 时 候 , 避 开 了 选择 
样本 受 损 区 域 。 在 物种 数量 一 定 的 情况 下 ， 随 着 每 
种 样本 数量 的 增加 , 识别 率 一 般 呈 上 升 趋势 , 并 稳 
定 在 茶 个 范围 ; 在 每 种 样本 数量 一 定 的 情 次 下 ,， 随 


不 同 物种 数量 时 每 种 识别 率 (% ) 


Identification rate of each species at diffierent number of species included 


3 spp. 4 spp. 5 spp. 6 spp. 7 spp. 8 spp. 
100 100 100 92.9 78.6 78.6 
75 75 91.7 75 75 66.7 
100 75 50 75 50 62.5 
一 62.5 50 62.5 37.5 25 
— — 57.1 57.1 14.3 28.6 
- - - 66.7 33.3 66.7 
- - - - 66.7 66.7 
- - - - - 0 
34 42 49 52 55 58 
91.2 81 75.5 75 56.4 55.2 


着 物种 数量 的 增加 , 识别 率 逐 渐 下 降 。 在 建 库 物 种 
总 数量 较 少 的 情况 下 ,每 种 样本 数量 较 多 时 ,物种 
的 识别 率 超过 90% , 甚至 达到 100% 。 而 在 建 库 物 
种 总 数量 较 多 的 情况 下 , 由 于 一 些 物 种 的 训练 集 样 
本 数量 太 少 , 无 法 满足 分 类 天 设 计 需 求 ， 导 致 测试 
结果 变 差 , 但 只 要 每 种 训练 集 样本 数量 大 于 5 幅 ， 
识别 率 均 大 于 80% 。 训 练 集 样本 数量 较 少 的 物种 
[如 有 (Z.) abbreviate, B. (B. ) carambolac, B. 
(2Z. ) caudata 和 B. (ZZ. ) scutellata 等 | 以 及 种 内 形 
态 差 中 大 的 物种 [| 如 B. (Z.) caudata, B. (8B.) 
dorsalis 和 B. (ZZ. ) scutellata 等 | 的 识别 率 较 差 。 

3 讨论 

本 文 的 研究 结果 表明 : 实 晶 科 果实 晶 属 昆虫 数 
字 图 像 自 动 识别 系 统 对 光照 不 均匀 、 样 本 变形 、 部 
分 破损 等 都 表现 出 了 良好 的 鲁 棱 性 , 日 识别 率 与 每 
种 训练 集 的 样本 数量 基本 正 相 关 , 与 建 库 物种 的 总 
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量 久 相关 。 以 上 结论 为 我 们 利用 “LBP 特征 + 
Adaboost 算法 ”进行 大 量 实 蝇 科 昆虫 的 物种 识别 提 
供 了 依据 , 即 尽 可 能 使 每 种 训练 集 的 样本 增加 到 合 
适 的 数量 。 此 外 , 本 系统 也 允许 使 用 者 将 已 经 正确 
命名 的 新 样本 加 入 到 训练 集中 , 从 而 进一步 提高 识 
别 率 , 因此 具 较 强 的 可 扩展 性 。 

但 目前 的 研究 没有 考虑 复杂 背景 去 除 的 问题 ， 
而 将 此 任务 交 给 用 户 自 行 处 理 , 这 虽然 提高 了 系统 
的 识别 性 能 , 但 却 降低 了 实际 应 用 效果 。 这 是 该 系 
统 未 来 需要 者 力 解 决 的 问题 之 一 。 此 外 ,虽然 本 系 
统 建立 在 已 有 模板 数据 库 基础 上 的 识别 结果 还 是 不 
错 的 , 但 建 库 标 本 量 较 小 物种 的 正确 识别 率 下 降 很 
快 。 因 此 在 下 一 步 工 作 中 , 我 们 将 尽 可 能 使 模板 数 
据 库 的 建 库 样本 量 增加 到 一 个 合适 的 数量 ,从 而 使 
本 系统 逐渐 达到 可 以 实际 应 用 的 目的 。 

本 系统 模板 数据 库 以 标本 图 像 为 主 , 然而 在 实 
际 生活 中 生态 图 像 也 有 有 自动 鉴别 的 巨大 和 需求。 生态 
图 像 更 为 复杂 , 需要 考虑 的 问题 更 多 。 在 进一步 的 
研究 中 , 作者 将 试 对 本 类 昆虫 生态 图 像 的 特征 及 其 
抽取 方法 进行 进一步 探讨 , 以 使 本 系统 具有 更 广 的 
应 用 范围 。 
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